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Abstract. In this paper, a meteoric plasma is analyzed 
from a physical viewpoint, with particular emphasis on 
its interaction with radio waves. The attention is drawn 
to some macroscopic characteristics of a meteoric plasma 
and it is shown that the electron-ion collision frequency 
is not negligible, as commonly thought. 
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1. Introduction 

When meteoroids enter the Earth's atmosphere, they cre- 
ate long and narrow trails of ionized gas, which can scatter 
radio waves. The meteor radio echo theory finds its roots 
in studies of the ionosphere made toward the end of the 



twenties (Skellet |1931 ). But only after the Second World 
War, and the development of military radar, the correla- 
tion between radio echoes and meteor trails became clear 
(Hey & Stewart 1946| ). The first experiments explicitly de- 
voted to meteor studies were carried out by Pierce (1947), 
who observed Draconids during the night of 9 to 10 Octo- 
ber 1946. 

The first theories on the interaction of radio waves with 
meteo rs were due t o Love ll & Clegg ( 194SD , Kaiser & Closs 
( |1952p , Herlofson ( |1951| ). Thereafter, between 1950 and 
1960, a lot of efforts in this field were undertaken, but after 
1960 the interest quickly decayed. The state of knowledge 



was well exposed in a classical paper by Sugar (1964) and 
the book by McKinley ( |1961| ). Toward the end of the eight- 
ies, the advent of digital technology renewed the interest 
in forward-scattering as a useful tool for communication 
channels over the horizon (see Weitzen & Ralston 1988| ). 
Today, meteor radars are widely used, even by a mateur 
astronomers, becaus e of lo w cost (Jenniskens et al. 1997 , 
Yrjola, & Jenniskens 1998 ). 

Actually, the two main theories for radio echoes are 
due to Poulter & Baggaley ( 1977]), w ho deal with back- 
scattering, and Jones & Jones (1990a, 1990b, 1991), who 



deal with forward-scattering. We can note that both the- 
ories do not consider the meteor as a plasma, but as a sim- 
ply ionized gas, with a negligible collision frequency. This 



is an assumption quite common in work on radio echoes 
from meteors, except for Herlofson ( 1951 ), who first con- 
sidered the meteor as a plasma. 

Concerning the difference between back-scattering and 
forward-scattering, it is worth noting that the two systems 
use different types of radio waves. A forward-scattering 
radar uses a continuous sine wave, while a back-scattering 
one uses a pulsed wave. This influences the mathemati- 
cal approach to the problem, but also physical theories. 
The pulse shape of the back-scattering radar can be rep- 
resented by a sum of several components with different 
frequencies. It follows that, during the propagation, the 
various components tend to change phase with respect to 
one another, which leads to a change in shape of the pulse. 
The dispersion relation thus depends on the frequency and 
can be expressed in a Taylor series. A detailed treatment 
of the pulse propagation in dielectrics is not the subject of 
this paper: the interested reader can find useful mathemat- 
ical tools in, for example, Oughstun & Sherman ( 1997 ). 
Here we want only to underline that we cannot simply 
consider the signal emitted by a back-scattering radar as 
sinusoidal: thus, the extension of the back-scattering echo 
theory to forward-scattering is not correct. 

The purpose of this paper is to settle some basic con- 
cepts in meteor physics and we leave a "cook-book" ap- 
proach to other articles. The interaction of sine waves (for 
the sake of simplicity) in the radio frequency range with 
meteoric plasma is investigated. We will see that the com- 
mon assumption about meteors as a collisionless ionized 
gas does not have any physical ground. 

2. Meteoric plasma 

As a meteoroid enters the Earth's atmosphere, it collides 
with air molecules. At the heights where most meteors 
ablate, the mean free path of the air molecules is about 
0.1 — 1 m. On the other hand, common meteoroid dimen- 
sions are of the order of 10 -3 — 10~ 2 m. This means that 
there is no hydrodynamic flow around the meteoroid and 
single air molecules impact on the body. If we consider a 
meteoroid a typical geocentric speed of 40 km/s, it can 
be found that air molecules impinge on the body with 
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the same speed. The kinetic energy is about 1.3 • 10~ 18 J 
(8 eV) per nucleon: a nitrogen molecule then has an en- 
ergy of about 3.7- 10~ 17 J or 230 eV. The impact energy is 
readily transformed into heat, which makes atoms evapo- 
rate from the meteoroid. The collisions between free atoms 
and air molecules produce heat, light and ionization, i.e. a 
meteor. Since this transformation occurs throughout the 
flight, the meteoroid atoms are dispersed in a cylindrical 
channel along the path. The electron line density is pro- 
portional to initial mass of the meteoroid, because the air 
mass involved is negligible when compared to the mete- 
oroid mass. 

After the escape, the first collisions of meteoroid atoms 
with air molecules take place at a distance of about one 
mean free path from the meteoroid path. It is useful to 
consider only the first collision to be important for ion- 
ization. This explains why the radio echo quickly rises to 
maximum amplitude and then slowly decays. At the mo- 
ment of creation, all electrons are thus located inside a 
cylinder with a radius of about one mean free path. 

It is possible to calculate the Debye length, a parame- 
ter that allows us to establish if the meteor is a plasma or 
simply an ionized gas. If the Debye length is small when 
compared with meteor characteristic dimensions, then it 
is possible to speak of plasma, i.e. a gas where the electro- 
static energy exceeds the thermal energy. In this case, if 
the thermal energy produces deviations from charge neu- 
trality, a strong electric field arises in order to restore the 
charge neutrality. On the other hand, if the meteor char- 
acteristic dimensions are small compared to the Debye 
length, this means that the thermal energy exceeds the 
electrostatic energy and there is no charge neutrality. In 
this case we have a simple ionized gas and we cannot speak 
of a plasma. It is important to underline the difference 
between a plasma and an ionized gas: a plasma has some 
macroscopic properties, such as the Langmuir frequency, 
which are absent in ionized gas. 

As known, the Debye length is obtained by equating 
the thermal energy to th e elec trostatic energy and equals 
(see Mitchner & Kruger [19731) : 



e a kT 



(1) 



where eo [F- m _1 ] is the vacuum dielectric constant, k [J- 
K _1 ] the Boltzmann constant, e [C] the elementary elec- 
tric charge, T [K] the temperature and n e [m -3 ] the elec- 
tron volume density. For the sake of simplicity, in this pa- 
per, the electron volume density is used, although in radar 
studies on meteors the electron line density is used. How- 
ever, it is easy to obtain the electron line density noting 
that the meteor trail is like a long circular cylinder. 

Common radio meteors are characterized by electron 
volume densities between 10 11 m -3 and about 10 20 m -3 



1995). By substituting these values in Eq. ([!]) it is possi- 
ble to calculate the maximum value of the Debye length, 
which is about 0.01 m. Comparing it to the lower value of 
the meteor initial radius (0.1 m), it is possible to say that 
the meteor is a plasma and not an ionized gas. 

Since any slight distortion of the plasma from a con- 
dition of electrical neutrality gives rise to strong restoring 
forces, we have to consider how fast these forces act. Writ- 
ing the equation of motion for the electrons, it is possible 
to find that they oscillate with a characteristic angular 
frequency (see Mitchner & Kruger 1973): 



e m e 



(2) 



where m e [kg] is the electron mass. It is possible to cal- 
culate Eq. (g) for any type of charged particle, but as 
electrons move faster than ions, they give the main con- 
tribution to plasma frequency and other contributions can 
be neglected. The characteristic angular frequency is often 
called the Langmuir frequency. 

Now, it is possible to understand the different be- 



(Sugar 1964) and temperat ures b etween 1000 and 5000 K 
(Bronshtcn 1983, Borovicka 1993, Borovicka & Zamorano 



haviour of radio echoes. Already Kaiser & Closs (1952) 
noted a different behaviour of meteors due to electron line 
density. They introduced the names underdense and over- 
dense to characterize meteors with electron line densities, 
respectively, below or above the value 2.4- 10 14 m _1 . In un- 
derdense meteors, the electron density is sufficiently weak 
to allow the incident wave to propagate along the ionized 
gas and the scattering is done independently by each in- 
dividual electron. On the other hand, when the meteor is 
overdense, the electron density is sufficient to reflect to- 
tally the incident wave. It has been usual to resort to a 
simple model in which the trail is regarded as a totally 
reflecting cylinder. 

Now, comparing the electromagnetic wave frequency 
to the plasma frequency (Fig. Q) we can distinguish two 
regions depending on whether the plasma frequency is 
higher or lower than the radar frequency. In the first case, 
the charges in the plasma have sufficient time to rear- 
range themselves so as to shield the interior of the plasma 
from the electromagnetic field (overdense meteor). In the 
second case, since the radar frequency is higher the than 
plasma frequency, the incident wave can propagate along 
the plasma (underdense meteor). Moreover, space charges 
can appear. 

Thus, a first method to distinguish underdense from 
overdense meteors, is by equating plasma frequency to 
radar frequency. We stress that this is valid for forward- 
scattering radar, which use continuous waves, while back- 
scattering radar see the meteors in a different way, because 
they use short pulses. 

3. General equations 

In a plasma the magnetic induction B and the magnetic 
field strenght H are almost always treated by the same re- 
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Fig. 1. Plasma frequency as a function of electron vol- 
ume density. Typical radar frequencies are also shown for 
comparison. SI units are used. 

lationship as in free space. In order to account for polariz- 
ability associated with bound electrons of neutral particles 
and ions, it is necessary to calculate the dielectric constant 
of the plasma. However it is possible to see that, in the 
radar frequency range, the contribution due to polariz- 
ability is negligible. Then, the relation between the elec- 
tric field strenght E and electric induction D can also be 
treated as in free space. The generalized Ohm law serves 
to relate J and E: 



J = aE = 



m e (v — iuj) 



(3) 



where v [s _1 ] is the mean collision frequency and a [S- 
m _1 ] is the electric conductivity. The meteoric plasma can 
be considered an isotropic medium since the electron gy- 
rofrequency is much less than the radio wave frequency 
and thus the effect of the geomagnetic field can be ne- 
glected. 

As stressed in Sect. 1, forward-scattering radars use 
continuous waves: taking into account that the time de- 
pendence of the electromagnetic field is e~ luJt , and using 
basic equations and standard identities, it is possible to 
combine Maxwell's equations in the following form: 

V 2 E + (uj 2 fi e + iuj^ a)E = V(V • E) (4) 

V 2 Jf + (w 2 ^ £o + iu/M)CT)H = V(er - iuoe ) x E (5) 

These are equations for the electric and magnetic fields 
in a meteoric plasma. It must be noted that Eqs. (^) and 
(||) are useful in the underdense region, where there are 



space charges. On the other hand, in the overdense region, 
when the plasma frequency is higher than the radar fre- 
quency, charge neutrality is present and then, the right 
side term of Eq. (|4|) vanishes, according to Gauss law, 
when the net charge density is zero: 



V 2 E + (w 2 // eo + iuiVo<r)E = 



(6) 



The solutions of Eq. @ are nonuniform harmonic 
plane waves of this type: 

E = E e l(k - r - UJt) (7) 

where the wave vector in Eq. (^) has the form: 

k = (3 + ia. (8) 

The two real vectors (3 and a. generally point in different 
directions. If we substitute Eqs. (0) and (||) in Eq. (||), it 
is possible to find that: 



k 2 = cj 2 /io£o + ioJUo® 



(9) 



Taking into account Eqs. (g) and and recalling that 
figeo = c~ 2 (c: light velocity in vacuum), it is possible to 
rearrange Eq. <M): 



(^) 2 (^) 



i+(^) 2 



(10) 



The quantity n — kr + ikj is often identified as the 
complex dielectric constant for the medium. From Eqs. (JsJ) 
and ([n]), when vectors (3 and a point to the same direc- 
tion, it is possible to define the real numbers a and (3 as 
the attenuation and phase constants respectively: 



a = — ., 

c V 2 



c 



KB. 



(11) 



(12) 



If the collision frequency is negligible, i.e. v << u>, then 
Eq. ( |l0| ) can be reduced to: 



{'-e» 2 } 



(13) 



that is, the equation used until now. If u > lo p , then k is 
a real number and the incident wave propagates without 
attenuation. On the other hand, if ui < lj p , then A: is a 
purely imaginary number and the incident wave is totally 
reflected. 
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Fig. 2. Typical radio echo from a meteor trail. The di- 
mensions are exaggerated. For an explanation of symbols, 
see text (Sect. 4). 

4. Collision frequency 

In all work on the theory of radio echoes from meteor 
trails, the collision frequency is considered 2 or 3 orders 
of magnitude lower than the radar frequency (a detailed 
analysis of past radio echo theories can be found in Fos- 



chini 1997). This is the typical collision frequency of elec- 



trons with air molecules at the heights where meteors ab- 
late. Work on collisions in meteor trails generally deals 
with ionization and excitation, in order to know processes 



during tr ail for mation (Massey & Sida 1955, Sida 1969 



Baggaley 198C). Other authors are interested in diffu- 
sion and thus study attachment, recombination and other 
chemical reaction s betw een the atmosphere an d me teoric 
plas ma (Ba ggaley 1972j , B aggaley & Cummack 1974 , Bag- 
galey 198C| , Jones & Jones 1990c ). In Fig. [2] a typical over- 
dense radio echo from a meteor trail is plotted. When the 
"flat top" (P-EP) is reduced to zero, we have an under- 
dense echo. Collision processes during the trail formation 
(from Start to Peak, S-P) are ionization and excitation, 
caused by air molecules impinging on the meteoroid. Dur- 
ing echo decay, from End Peak (EP) to End Decay (ED), 
attachment and recombination are dominant processes, 
owing to the diffusion of the plasma in the surrounding at- 
mosphere. During the plateau (P-EP) the trail can be con- 
sidered in local thermodynamic equilibrium, when matter 
is in equilibrium with itself, but not with photons (Mitch- 
ner & Kruger 1973). The plasma in the trail is transpar- 
ent to radiation at some optical frequencies, which escapes 
from the meteor and form the light we observe. 

For a plasma in local thermodynamic equilibrium it is 
still possible to employ the Boltzmann and Saha equations 



(sec Mitchncr & Kruger 1973), even though complete ther- 
modynamic equilibrium does not prevail. Now, thermal 
ionization becomes the main process and metals, with low 
ionization energy, drive this phase. It must be noted that 
in meteoroids there are some few per cent of alkaline and 
alkaline-earth metals, such as Na, K, Ca and Mg. These 
data are obtained from studies on meteorites (Mason 1971 , 



Millman 1976, Wasson 1985) and on bright fireball spectra 
(Bronshtcn 1983, Borovicka 1993] , Borovicka & Zamorano 
|1995 ). Because of their low ionization energy, at typical 
radar meteor temperatures (1000-5000 K), these metals 
easily ionize and thus, they are the main contributors to 
thermal ionization. A small percentage of one of these met- 
als suffices to produce a huge amount of electrons. 

In the calculation of collision frequency of electrons 
with other particles, it must be taken into account that, 
because of the electrostatic field, ions have a cross section 
larger than that of atoms and molecules. Thus, electron- 
ion collisions must be considered, and not electron-air 
molecule collisions, since the former are more frequent 
than any other. 

Moreover, in calculating the electron-ion collision fre- 
quency we have considered potassium, because it has a 
low ionization energy (4.34 eV). This choice could be 
questionable because potassium is scarcely present in me- 



teor spectra (Bronshten 1983) or even absent (Borovicka 



1993, Borovicka & Zamorano 1995). However, this metal 



is present in almost all metcoritic specimens (Mason 1971) 
and other studies on meteor shower have reported its 



presence. Specifically Goldberg & Aikin (1973), using a 
rocket-borne ion mass spectrometer, detected the pres- 
ence of potassium ion in /3-Taurids. The absence of potas- 
sium in meteor spectra is due to the low efficiency of light 
emission processes in this type of atom: it is well known 
that the most intense spectral line of potassium is in the 
infrared range (766.49 nm). It is worth noting that also 
iron and sodium, two of the commonest elements in me- 
teor spectra, could be absent in some cases (see Millman 
1976). Therefore, we can consider a potassium percentage 
of about 1% of meteoroid mass. 

It is very difficult to find experimental values for 
electron-ion collision cross sections of various species, in 
the temperature range of meteors. Experimental studies 
are mainly devoted to excitation and ionization cross sec- 
tions, in order to compare results with data from spectra 
(Ncff |1964[ Bo itnott & Sava ge |1970| , |l97l| , |1972| , Savage & 
Boitnott 1971). Rosa (1987) gives collision cross sections 
of K, K + and some atmospheric gases, such as O2, in the 
temperature range from 2000 K to 3500 K. It is important 
to note that the K + cross section is about 3 (three) orders 
of magnitude larger than those of other species. 

We have further supposed that only a single ionization 
is possible, i.e. we have a reaction such as: 



K K+ + e" 



(14) 



We would like to stress that we are only interested in 
the analysis of meteoric plasma in a steady state condi- 
tion (during the plateau) and we do not actually consider 
echo decay (recombination and other chemical processes) 
or trail formation (collisional ionization and excitation). 




Fig. 3. Electron-potassium ions collision frequency as a 
function of electron density. Temperature 3000 K. Typical 
radar frequencies are also shown for comparison. 

Thus, it is possible to use a simplified model of collision 



frequency (Mitchner & Kruger 1973), viz 



(15) 



where Q e i is the electron-ion collision cross section and 
v e i is the electron mean velocity with respect to ions (it is 
assumed to be about equal to the electron mean thermal 
velocity), and rij is the ion volume density. In our case, 
owing to Eq. (fbi]), rx, = n e . 

Now, it is possible to observe that v e i cannot be negligi- 
ble anymore, especially for a high electron volume density 
(Fig. H). If we consider the contributions of other species, 
it is necessary to sum the various collision frequencies ob- 
tained using Eq. ([l5|). 

Now, in the overdense regime, we have to consider the 
dispersion relation in Eq. ( |To|) and not in Eq. (|l3|). It is 
worth noting that Eq. (13) is still valid in the underdense 
regime because at high electron densities only the colli- 
sions are not negligible. 

5. Forward scatter of radio waves from meteoric 
plasma 

It is possible to define, from Eq. ([l^), a critical electron 
density in order to separate the underdense from the over- 
dense regime and this occurs when: 



1 



1 + (S) 2 







(16) 



In order to make some example, we can assume a 
mean temperature of 3000 K. If we consider the CNR 



Fig. 4. Wave vector modulus as a function of electron 
density above critical value. Temperature 3000 K; Radar 
frequency 42.7 MHz. Line (A): with collisions; line (B): 
without collisions. 



radar facility (Cevolani et al. 1995), which has a radar fre- 
quency of 42.7 MHz, we obtain a critical electron density 
of n e = 2.3 • 10 13 m~ 3 . On the other hand, if we consider 
the M. de Meyere's radar (see C. Steyaert, Radio Meteor 
Obs. Bulletin, ftp://charlic.luc.ac.be/pub/icaros/rmob/), 
which has a frequency of 66.51 MHz, we obtain a criti- 
cal electron density of n e = 5.5 • 10 13 m~ 3 . Below this 
critical density (underdense meteors), the wave number is 
real and the incident wave propagates into the meteoric 
plasma with negligible attenuation. We obtain well-known 
results about underdense meteors. Above this value, waves 
which were previously excluded (see Eq. (|l3|)), can now 
propagate, but are strongly attenuated (Figs. f| and ||). 
From Fig. ||, it is possible to see that the presence of a 
collision frequency comparable with the radar frequency 
determines a rise of the real part of the wave vector modu- 
lus. Physically, this behaviour may be understood on the 
basis that collisions subtract energy from plasma oscil- 
lations, allowing the incident wave to penetrate, even if 
strongly attenuated. 

We can observe two types and two sub-types of be- 
haviour. We have two main classes of meteors, overdense 
and underdense, according to whether the plasma fre- 
quency is higher or lower than the radar frequency respec- 
tively, depending on the initial electron volume density. 
Moreover, overdense meteors are divided into two sub- 
classes (see Fig. |^): for an electron volume density between 
the critical density and about 10 17 m~ 3 , the real part of 
the wave vector is negligible and the incident wave is to- 
tally reflected (overdense I). For a higher electron volume 
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Fig. 5. Attenuation and phase constants, taking into ac- 
count collision frequency, as a function of electron den- 
sity above critical value. Temperature 3000 K; Radar fre- 
quency 42.7 MHz. 

density, collisions allow the propagation of the incident 
wave, even if with a strong attenuation (overdense II). 

Total reflection, in the overdense I range, allow us to 
make some useful approximations in calculating the at- 
tenuation of radio waves after the forward-scattering. If 
the plasma had a definite boundary and a uniform elec- 
tron density, then reflection at its surface would be simple 
because the gas would act as a dielectric with a complex 
dielectric constant. But a meteor trail does not have a def- 
inite boundary and thus, the incident wave penetrates a 
little into the plasma before reaching the density neces- 
sary to allow total reflection. Reflection occurs gradually, 
as in a mirage. 

Taking into account this fact, it is possible to make 
some approximations. We can consider a simple geometry, 
as shown in Fig. ^, and then use the definition of the 
attenuation a in decibel units: 



a - 10 - log jfF [dB] 



(17) 



where subscripts i and r stand for incident and reflected 
wave. We substitute Eqs. (0) and (§) in Eq. @ and, tak- 
ing into account that the amplitude of a totally reflected 
wave is equal to the amplitude of the incident wave, we 
can obtain an attenuation value of about a — — 20aZloge, 
where I is the path of the wave into the plasma. From 
Fig.@ we can see that: 
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Fig. 6. Forward-scatter reflection geometry. 

where 6 is the penetration depth and <p is the incidence 
angle. Now, if we consider something similar to the "skin 
effect" in metals, we have 8 = 1/a. Then, Eq. (0) be- 
comes: 



-401oge -17.36 



47.36 sec [dB] 



(19) 



The attenuation is simply a function of the angle of 
incidence and this is compatible with results obtained by 



Forsyth & Vogan (1955), who predicted an attenuation 
proportional to sec 4>. It is necessary to stress that Eq. ( 19 ) 
is valid only for type I overdense meteors, where total 
reflection is allowed. 

In order to make some comparison with experimental 
data, we can consider a sample of radio echoes obtained 
from the CNR forward-scattering radar facility (e.g. Fos- 

1995| ). Attenuation of 



chini ct al. 1995 



Porubcan et al. 
reflected waves, in the overdense regime, ranges from -67 
to -47 dB. By using the radar equation (Kingsley & Que- 



gan 1992), attenuation values range from -155 to -125 dB, 
which gives a strong discrepancies with what is observed. 
To obtain from Eq. ((l|) the attenuation values requested, 
it is necessary to use incidence angles from 68° to 75°. This 
is a very good approximation because the CNR radar has 
the main beam with about 15° elevation angle: then 75° 
is just the complementary angle. 

For overdense type I meteors, i.e. those with a plateau, 
Eq. ( |l9[ ) can be used to calculate the meteor height. Tak- 
ing into account the CNR radar geometry (transmitter- 
receiver distance is about 700 km, see Cevolani et al. 
|1995 ), it is possible to obtain a range of overdense type 
I meteor heights from about 94 to 141 km. This is very 
important, because the height is obtained without tak- 
ing into account the diffusion coefficient, which is very 
uncertain. This requires further investigations, which are 
currently being carried out. 

6. Conclusions 



In this paper, we settled some basic concepts in meteor 
fity physics. We have dealt with the interaction of sine waves, 
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in the radio frequency range, with meteoric plasma. At- 
tention is drawn to some macroscopic characteristics of a 
meteoric plasma and it is shown that the electron-ion col- 
lision frequency is not negligible, as commonly thought. 
It is possible to define two meteor classes (ovcrdensc and 
underdense) according to whether the plasma frequency 
is higher or lower than the radar frequency respectively. 
Overdense meteors are divided into two sub-classes (I and 
II) , depending on the ratio between collision frequency and 
radar frequency. 

Taking into account that the meteoric plasma does 
not have a definite boundary, a simple formula for the 
calculations of radio-wave attenuation after the forward- 
scattering is also presented. This formula allows us to cal- 
culate the meteor height, without taking into account the 
diffusion coefficient. 

Further questions can be put: in this work, potassium 
ion is considered, but other studies can be carried out in 
order to know the impact of other alkaline and alkaline- 
earth metals, such as Na, Ca and Mg. 

Acknowledgements. Author wishes to thank G. Cevolani, 
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